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ABSTRACT 

Careful monitoring of harmonically bound (or as a limiting case, free) 
masses is the basis of current and future gravitational wave detectors, and of 
nanomechanical devices designed to access the quantum regime. We analyze 
the effects of stochastic localization models for state vector reduction, and of 
related models for environmental decoherence, on such systems, focusing our 
analysis on the non-dissipative forced harmonic oscillator, and its free mass 
limit. We derive an explicit formula for the time evolution of the expectation 
of a general operator in the presence of stochastic reduction or environmen- 
tally induced decoherence, for both the non-dissipative harmonic oscillator 
and the free mass. In the case of the oscillator, we also give a formula for 
the time evolution of the matrix element of the stochastic expectation den- 
sity matrix between general coherent states. We show that the stochastic 
expectation of the variance of a Hermitian operator in any unraveling of the 
stochastic process is bounded by the variance computed from the stochas- 
tic expectation of the density matrix, and we develop a formal perturbation 
theory for calculating expectation values of operators within any unraveling. 
Applying our results to current gravitational wave interferometer detectors 
and nanomechanical systems, we conclude that the deviations from quan- 
tum mechanics predicted by the continuous spontaneous localization (CSL) 
model of state vector reduction are at least five orders of magnitude below 
the relevant standard quantum limits for these experiments. The proposed 
LISA gravitational wave detector will be two orders of magnitude away from 
the capability of observing an effect. 
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1. Introduction 



Testing whether quantum mechanics is an exactly correct theory, or is an approxi- 
mate theory from which there are small deviations, is a subject of current theoretical and 
experimental interest. Significant bounds have been set [1] on deterministic, nonlinear mod- 
ifications of the Schrodinger equation [2], and such modifications are also theoretically disfa- 
vored because they have been shown [3] to lead to the possibility of superluminal communi- 
cation. On the other hand, stochastic modifications to the Schrodinger equation have been 
extensively studied [4] as a way of resolving the measurement problem in quantum mechan- 
ics, and are known to be theoretically viable. This raises the question of what bounds on 
the stochasticity parameters are set by current experiments, and what degree of refinement 
of current experiments will be needed to confront, and thus verify or falisfy, the stochastic 
models. 

The most extensively studied stochastic models are those based on the concept of 
localization [4,5], in which a stochastic, Brownian motion couples to the system center of 
mass degree of freedom. Weak bounds on the stochasticity parameters for this type of model 
can already be set [6] from experiments [7] observing fullerene diffraction, and stronger (but 
far from definitive) bounds will be set [8] by a recently proposed experiment [9] that aims to 
coherently superimpose spatially displaced states of a small mirror attached to a cantilever. 
Our aim in this paper is to analyze the effects of stochastic localization on another class of 
precision experiments, involving the careful monitoring of massive objects in gravitational 
wave detectors, and of microscopic oscillating beams in nanomechanical experiments. To 
this end, we analyze the stochastic Schrodinger equation for a non-dissipative forced har- 
monic oscillator, focusing particular attention on the effects of the stochasticity terms on 
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the quantum non-demolition variables of the oscillator. We also derive analogous formulas 
for the limiting case of a free mass, correcting a factor of 2 error in previous formulas given 
in the CSL literature. Because the stochastic expectation of the density matrix in the local- 
ization model obeys a differential equation used as a model for environmental decoherence, 
our results are also relevant to the study of decoherence effects on both the forced oscillator 
and free mass systems. Analyzing various experiments using our results, we conclude that 
for the parameters of current gravitational wave detectors and nanomechanical beams, only 
weak bounds will be set on the CSL model stochasticity parameters. The proposed LISA 
gravitational wave detector should do better, but is still not expected to see an effect. 

This paper is organized as follows. In Sec. 2 we give the basic stochastic Schrodinger 
equation to be analyzed, the corresponding pure state density matrix equation, and the 
simpler equation for the stochastic expectation of the density matrix (which is the usual 
mixed state density matrix). The latter equation, we note, is also used as a model for 
environmental decoherence effects, and so its solution is of particular interest. We also review 
briefly the basic ideas of quantum non-demolition measurements, leading to the identification 
of the non-demolition variables of the forced harmonic oscillator. In Sec. 3 we give results 
for the time evolution of expectations of the non-demolition and other low order variables 
of the forced oscillator. For comparison with the zero frequency limit of the oscillator, we 
give in Sec. 4 analogous results for a free mass, rederiving (and correcting) results already 
in the literature. In Sec. 5 we give formulas for the time evolution of stochastic expectations 
of general operators for the forced non-dissipative oscillator and for its free mass limit, 
and additionally derive a formula for transition amplitudes of the oscillator, giving results 
that also apply to environmental decoherence effects. In Sec. 6 we consider stochastic 
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fluctuations, and show that expectations of variances of observables can be bounded using 
our earlier calculations proceeding from the expectation of the density matrix. In Sec. 7 
we set up a formal perturbative procedure for calculating stochastic fluctuation effects, and 
use the leading order results to interpret the inequality derived in Sec. 6. Finally, in Sec. 
8, we apply our results to make estimates for the effects of CSL models in gravitational 
wave detection and nanomechanical resonator experiments. In Appendix A we review some 
Ito calculus formulas, and in Appendix B we relate the formalism used in the text to the 
Lindblad density matrix evolution equation. 



We start our analysis by considering a massive harmonic oscillator in one dimension, 
which in the three-dimensional case will describe the dynamics of one center-of-mass degree 
of freedom. The oscillator Hamiltonian is taken as 



with uj the oscillator angular frequency, d(t) a complex c-number driving term, and a, a* 
annihilation and creation operators obeying [a, a*] = 1. These operators are related to the 
oscillator mass m, coordinate q, and momentum p, by 



2. Basic formalism: one dimensional oscillator 



H = %u(a ] a + -) + d(t)a) + d(t)a 



(1) 



a ={mui/2h) 2 (q + ip/mu) 




(2a) 



q =<j(a + a') 



a = {Ti/2mui) 



and the number of quanta iV in the oscillator is given by 



N = a f a 



(26) 
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Discussions of quantum non-demolition experiments involving oscillators [10] also introduce 
the quantities 

X 1 —q cos out — (p/muj) sin cut , 

(3a) 

X 2 =q sin cut + (p/ moj) cos out , 

from which one easily finds 

q + ip/mu =(X ± + iX 2 )e- iujt , 
q- ip/mu =(Xi -iX 2 )e iwt , 

(36) 

X 1 =a(ae luJt + oV*"*) , 

X 2 = - ia(ae luJt - ah~ iwt ) . 
Hence X^ 2 are quantum mechanical analogs of the classical amplitude of the oscillator, and 
when the external driving term d(t) is zero they are conserved, as is the occupation number 
N. Because these quantities are constants of the motion in the absence of external forces, 
measurements of them, while introducing uncertainties into the conjugate variables (which 
are the phase in the case of N, X 2 in the case of X i: and X 1 in the case oiX 2 ), do not feed 
the uncertainties in the conjugate variables back into the time evolution of the measured 
variable. Hence the variables N, Xi, and X 2 can in principle be measured to an accuracy 
not limited by the uncertainty principle, and are called "quantum non-demolition" variables. 

Letting \ip t ) be the oscillator wave function at time t, the standard Schrodinger 
equation is 

d\i/> t ) = -(i/h)Hdt\^ t ) . (4a) 

We shall be interested in this paper in a class of models [11] for state vector reduction, which 
modify Eq. (4a) by adding stochastic terms to the Schrodinger equation. Specifically, we 
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shall consider the evolution equation 



d\A) 



- l -Hdt + ^]{q- (q)) dW t ~\{q- (q)) 2 dt 



, (46) 



where H is given by Eq. (1), and (q) = (ipt\ ( l\i J t) is the quantum mechanical expectation of 
the position operator q of the oscillator. Introducing the pure state density matrix p(t) = 
\i])t) (i]>t\i we can a l so write (q) = Trgp(t). The stochastic dynamics is governed by a standard 
Wiener processes W t , defined on a probability space (VL,F, P). Using the rules of the Ito 
calculus (see Appendix A), the density matrix evolution corresponding to Eq. (4b) is 

i 1 

dp= -j\H,f>]dt- -rj[q, [q,p}}dt + y/fj\p,\p,q]]dW t ■ (5a) 

Since this evolution equation obeys {dp, p} + (dp) 2 = dp, it preserves the pure state condition 
p 2 = p. When statistics are accumulated by averaging many runs of an experiment, the 
relevant density matrix in the stochastic case is the ensemble expectation p = E[p], giving 
the mixed state density matrix which obeys the ordinary differential equation 

dt h 2 (5ft) 
i i 

= - j\H,p\- -r]a 2 [a + a\[a + a\p\] . 
This equation is of particular interest because (with a different value of the parameter 77) 
it is also used [12] as a simple model for environmental decoherence effects. The calcula- 
tions of this paper focus on analyzing Eqs. (5a) and (5b) for the forced harmonic oscillator 
Hamiltonian of Eq. (1). 

3. Stochastic expectations of oscillator observables 

We begin by considering the evolution equation of Eq. (5b) for the mixed state 
density matrix p. Letting B be any time-independent operator, let us denote by ((B)) the 



expectation computed with the mixed state density matrix, 

((B)) =Tr pB . (6a) 

For the time evolution of this expectation, we then find 

d((B)) =Ti dp B 
dt dt 

=T±B(-±[H,p]-± V [q,[q,p]]) (66) 
=Tr( - l -[B, H] - ^ W 2 [[B, a + a^a + a^p , 
where we have made repeated use of cyclic permutation under the trace. Let us now make 
successively the choices B = a,a^ , aa, aW , a^a, aa^ = 1 + a^a, corresponding to all quantities 
linear and quadratic in the creation and annihilation operators. Then evaluating the single 
and double commutators in the final line of Eq. (6b), a simple calculation gives for the two 
linear operators, 



d i 

— Trpa = — iuTrpa — —d(t) , 

UL ft 

d i — 

— Tr pa f =iwTrpa! + -d(t) , 

UL f V 

and for the four quadratic operators 

d % , , « q 

— Tvpaa = — 2iuoTipaa — 2—d(t)Tvpa — r\a 

UL I V 

d %- 

— TYpaV =2^TrpaV + 2-d(t)Trpa f - r]a 2 

UL I L 

—Tipa)a = — ^-d(t)Trpa^ + — d(t)Trpa + r/a 2 

UL I i lb 

d t d + 

— Trpaa 1 =— Tr pa' a . 

These equations can be immediately integrated to give 

Trp(t)a =e~ iujt [Trp(0)a - % - f dud(u)e^ 
L n Jo 

Trp(t)a f =e tujt [Trp(0)a f + % - [ dud{u) e - 1 ' 

n Jo 



(7a) 



(76) 



(8a) 



for the linear operators, and 



Tvp(t)aa = e - 2iwt Tvp(0)aa - dve 2iuv (2^d(v)Trp(v)a + r]a 2 ^j 

Trp(t)aV =e 2lu}t [Trp(0)aV + f dve~ 2wv (2^d(v)Trp(v)a^ - V a 
L Jo ^ n> 



(8b) 



Tip(t)a ] a =Trp(0)a t a - % - [ dv \d(v)Trp(v)<J - d(v)Trp(v)a\ + r]a 2 t , 

* Jo L J 

Trp(t)aa f =1 + Trp(t)a f a 

for the quadratic operators. We see that by substituting Eq. (8a) into Eq. (8b), we can 
reduce the expressions for the quadratic operators to quadratures. (Proceeding in a similar 
fashion, it is easy to see that given any polynomial P(a, aJ) of finite degree in the creation 
and annihilation operators, the expectation Trp(t)P can be reduced to quadratures; for an 
explicit formula constructed by generating function methods, see Sec. 5.) 

Rather than exhibiting the full expressions for the expectations of the quadratic 
operators, we note that what we are most interested in is calculating the change in these 
quantities, denoted by 5, arising from the "decoherence" term with coefficient 77 in Eq. (5b). 
From the fact that Eq. (8a) contains no terms proportional to 77, we see that there are no 
stochastic (or decoherence) effects on the linear operators, 

5Trp(t)a =0 , 

(9a) 

5Tip(t)a^ =0 , 

while the effect of the 77 term in Eq. (5b) on the quadratic operators is simply given by 

Tier 2 

5Tr p(t)aa = e~ lLOt sin out , 

LO 

2 

6Trp(t)aW = - ^-e^smcot , ( 96 ) 

LO 

5Tip{t)a ] a =5Trp(t)aa^ = rjaH . 
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Using the definitions of X lj2 given in Eq. (3a,b), we correspondingly find that 

5Trp(t)X 1 =5Trp(t)X 2 = , 
5Trp(t)X 2 =2r}a 4 (t- 



sin2a;t\ 
2u ) 
sin 2ut 



5Trp(t)X 2 =2 W ^t+ S -^) , (10) 
4na 4 

5Trp(t) (X ± X 2 + X 2 X 1 ) = — sin 2 ut , 

5Trp(t)[X 1 ,X 2 ]=0 . 

We note that these formulas are exact (not just approximations to first order in rj), since for 
all the operators B considered above, we have 

Tr p(t)B = Tr p(t)B\ v=0 + 5Trp(t)B . (11) 



4. The free mass limit 

According to Eq. (9b), the oscillator occupation number iV = a^a contains a term 
that grows linearly in time as r\o 2 t. Since the occupation number contribution to the oscil- 
lator energy of Eq. (1) is huN, and since a 2 = h/2mu from Eq. (2a), the oscillator energy 
contains a term that grows linearly in time as 

5E = 5Tip(t)H = ^ . (12) 

Because this formula is independent of the oscillator frequency u>, it should also correspond 
to the energy increase of an unbound mass m arising from the ?7 term in Eq. (5b). This can 
be calculated directly as follows. For an unbound mass in one dimension, the Hamiltonian 
is H = p 2 /2m, and the density matrix evolution is given by the first line of Eq. (5b). So we 
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have, by the same reasoning that led to Eq. (6b), 

^-STrpH = - Ti^-r][[p 2 , q], q]p 

at Am Mn x 

i (13a) 

= — Tr r][—2ihp, q]p = r/h 2 /2m , 

Am 

giving 

5TrpH = r]h 2 t/2m , (136) 

in agreement with the result calculated for the oscillator. This result is a factor of two 

larger than the one quoted in the CSL literature [13]; for instance, Taylor expansion of 

Eq. (3.36) of Ghirardi, Pearle, and Rimini (GPR) shows that for a uniform cube, their 7^ 

is the same as the parameter 77 used here, and so their formula of Eq. (3.38c), which states 

that ^((-Pj 2 )) = \l5ih 2 would correspond to ^((if)) = rjh 2 /Am, in disagreement with our 

result of Eq. (13a) and with the oscillator calculation of the preceding section. This error 

propagates through to Eqs. (3.41a) through (3.41c) of GPR, all of which are a factor of 2 

too small. Thus, in our notation, their results should read 

5Trpp 2 =r]h 2 t , 

STrp(pq + qp) =rjh 2 t 2 /m , (13c) 

5Trpq 2 =r]h 2 t 3 /(3m 2 ) . 
A rederivation of the second and third lines of Eq. (13c) will be given in the next section. 
[These equations were first given, with a different identification of the proportionality con- 
stant 77, in the GRW model [11]. Philip Pearle has rechecked the calculations in the paper 
of GPR, and finds that a factor of 2 error was made in going from their Eq. (3.36) to their 
Eq. (3.38c); when corrected, their equations agree with our results of Eq. (13c) above.] 
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5. Exact general formulas for decoherence effects on a forced 
non-dissipative harmonic oscillator and on a free mass 



We have seen in Eqs. (8a) and (8b) that the double expectations (stochastic and 
quantum, as defined in Eq. (6a)) of low order polynomials in the oscillator creation and 
annihilation operators can be reduced to quadratures. To show that this is a general result, 
let us consider the generating function 

K ap (t) =Tr(exp(«a t e^*) exp((3ae iujt )p(t)) 

00 00 n am (14a) 

^ ^ n\ ml w 

n=0 m=0 

from which one can extract the expectations of arbitrary normal ordered operators formed 
from a and a 1 " . To proceed, we shall need the generalization of Eq. (6b) to the case when 
the operator B has an explicit time dependence, which reads 

at v at ot 

/ rt¥i i 1 \ 

=Tr VM ~ h [B ' H] " ^ 2[[s ' a + ° t] ' a + a]] ) p ■ 

Applying this formula to Eq. (14a), with B = exp(aa^ e~' luJt ) exp(/3ae^*), the explicit time 
derivative on the right cancels the commutator term involving the free Hamiltonian Tiuoa^a 
(this is why we included an explicit time dependence in the definition of the generating 
function), leaving the simple differential equation 



Uae-'"'d(t) - /Je"'d(t)) - ir)<T 2 (ae-<"' - I3e'"'f 

it Zi 



K aP {t) . (15) 



Defining 



D(t) = I due iuju d{u) , D{t) = [ due- lu}U d{u) , (16a) 
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the integral of Eq. (15) takes the form 




2 



{a 2 e~ iwt + (3 2 e lu}t ) sinut + UaD(t) - f3D(t)) K aP (0) 



K a p(t) = exp a(3rja 2 t 



(166) 



with 



K aP (0) = Tr(e Qa VV(0)) 



(16c) 



Expanding this equation through second order in a and (3, one can verify that it agrees with 
the formulas of Eqs. (8a) and (8b), and so we have obtained the generalization of these ex- 
pressions to arbitrary normal ordered monomials in the creation and annihilation operators. 
Thus expectations of operators with respect to the density matrix of the decoherent forced 
oscillator can be explicitly calculated in closed form. As an example of particular interest, 
we note that the 77-dependent terms with the dominant time dependence for large times 
can be read off from the power series expansion of the first factor on the right hand side of 



Thus, the leading 77 dependence in Tr p{t)a) a at large times is i]a 2 t, in agreement with 
Eq. (8b), while the leading 77 dependence in Tip(t)a) a) aa is 2rfaH 2 . We will apply these 
results below to a discussion of the variance of N at large times. 

The same strategy that we have just followed can be used to find a generating function 
for the expectations of general polynomials in the operators q and p in the free particle case. 
Here the Hamiltonian is H = p 2 /(2m), and the equation to be solved is 



Eq. (16b), 




(16d) 



d((B)) 




dt 



(17a) 
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We consider now the generating function defined by 

K f a/3 = TrBp(t) = Tr[exp(a(q-tp/m))exp((3p)p(t)} . (17b) 
Using the fact that 

exp (— (i/Ti)tp 2 /(2m))gexp ({i/Ti)tp 2 /(2m)) = q — tp/m , (17c) 
we see that the terms dB/dt and — (i/h)[B, H] in Eq. (17a) cancel, so that we are left with 

±K^(t) = Tr(-^r ] [[B,q},q}p(t)) . (18a) 

Using now the identity 

exp (a(q — tp/m)) exp({3p) = exp(ccg) exp (p(/3 — at/m)) exp [c?%htj (2m)) , (186) 

the right hand side of Eq. (18a) is easily evaluated to give 

± V h 2 ((3-at/m) 2 K^(t) . (18c) 

Equations (18a) and (18c) now give a differential equation that can be immediately inte- 
grated, giving a result analogous in form to Eq. (16b), 

1 



K f ap {t) = exp 



Q r]h 2 t(3(3 2 - 3(3at/m + a 2 t 2 /m 2 ) 



with 

^i(0)=Tr(exp(ag)exp(^)p(0)) . (196) 

This equation gives a generating function from which the results of Appendix E of Ghirardi, 
Rimini, and Weber [11] and their extensions to higher order polynomials, can be readily 
extracted. In particular, expanding Eq. (19a) through second order in a and (3, one gets 
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for the leading 77 dependence of the expectations of quadratic polynomials in p and q the 
expressions given above in Eq. (13c). 

Returning to the harmonic oscillator, the same methods can be applied to the gen- 
erating function for general matrix elements of p(t), although the results in this case are not 
so simple. Let us define the generating function 

L af} (t) = Tr(exp(cmV iw *)|0)(0|exp(/3ae^*)p(i)) , (20a) 

where |0) is the oscillator ground state obeying a\0) = (0|a^ = 0. (With the inclusion 
of a normalization factor exp ( — |(|o;| 2 + |/9| 2 )), this expression gives directly the matrix 
element of p(t) between coherent states of the oscillator parameterized by a and (3.) When 
we take the time derivative of this expression, and apply Eq. (14b), we now find that there 
are additional terms where an a) multiplies |0) from the left, or an a multiplies (0| from the 
right. These can be converted to derivatives of L a/3 with respect to the parameters a and (3, 
and so we end up with the differential equation 



d r . . f i 

Jt L ^ = \n 



1 2 
'2" 



«-|) e ""-("-^) e " 2 



(206) 
L a p(t) , 



which corresponds to making the substitutions a — > a — (3 — * (3 — ^ in Eq. (15). Since 
the operators — a and ^ — (3 commute with one another, this equation can be formally 
integrated without requiring a time ordered product. Using 

— - p\ = e af3 d c~ ap 
15 



the result can be compactly written as 

L a/ 3 = e af3 exp 



d d 1 2 sincjt / _ iu!t ( d \ 2 iu!t ( d x 2 
rja t—— -r)o e 77^ +e 



dp da 2 iv \ \ dp J \da 



(216) 

with L° a p{t) the generating function in the absence of decoherence (that is, with 77 = 0), 
which is given by 

L%(t) = exp (-±\D(t)\ 2 + l -\D(t)a - D(t)P}^j Tre a ^ a+ ^ D ^ \0) (0\e a ^ h ^ p(0) . 

(21c) 

An alternative form of this result is obtained by introducing the Fourier transform of e~ ct/3 L° /3 
with respect to a and P, 

e~^L%{t) = j ' dp a dp p F(p a ,p p ,t)e iap " +i ^ , (22a) 

in terms of which Eq. (21b) takes the form 

L a p{t) =e al3 [ dp a d Pf3 exp - W 2 tPaPp + ^a 2 ^^ (e"^ + e lu)t p 2 a ) 

J L 2 " J (226) 

xF(p a , Pf3 ,t)e iap ° +i ^ . 
Thus, matrix elements of the density matrix for the decoherent forced oscillator can be 

explicitly (if formally) expressed in terms of matrix elements of the oscillator in the absence 

of decoherence. 

We have seen that exact results can be obtained for a number of properties of the 
density matrix evolution equation of Eq. (5b). This might have been suspected from the fact 
that earlier work [8] has shown that this equation leads to an exactly solvable expression for 
the fringe visibility in a proposed mirror superposition experiment described by an oscilla- 
tor Hamiltonian. More general density matrix evolution equations for a damped harmonic 
oscillator have been discussed in the literature [14]. When additional decoherence terms 
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of the form [c\q + C2P, [c\q + C2P, p}} (for general constants 01,2) are added to the density 
matrix evolution equation for the forced oscillator, an explicit result for K a p{t) generaliz- 
ing Eq. (16b) can still be easily obtained. When dissipative terms proportional to a linear 
combination of i[q, {p, p}} and i\p, {q, p}}, with { , } the anticommutator, are added to the 
density matrix evolution equation, the differential equation for K a p(t) contains terms involv- 
ing d/da and d/d{3, and we then can no longer obtain an explicit formula for the expectation 
of the generating function analogous to Eq. (16b). Such dissipative terms are included in 
the evolution equations discussed in refs. [14], where some exact results are obtained. We 
remark, however, that for mechanical or electrical systems with a very high quality factor 
Q, it can be a useful first approximation to neglect classical damping in studying stochastic 
reduction and decoherence effects, as done in the analysis of this paper. [For the benefit of 
the reader familiar with the Lindblad form of the density matrix evolution equation, we give 
in Appendix B its relation to the commutator/anticommutator structures discussed here.] 

6. Bounds on variances for unravelings 

So far we have studied quantum expectations of physical quantities in the mixed state 
density matrix p obtained as the stochastic expectation of the pure state density matrix 
p that obeys Eq. (5a). In any given run of the physical process (or "unraveling" in the 
stochastics literature parlance), the quantum expectation of a physical quantity represented 
by a non-stochastic operator B will be governed by TrpB. As before, let us use the notation 
(■ • •) to denote expectations formed with respect to p, and the notation ((• • •)) to denote 
expectations formed with respect to p = E[p]. Then by linearity we evidently have 

((B)) = E[(B)] . (23) 
17 



We shall now show that the variances corresponding to the single and double averages are 
related by an inequality. Let 

((AB) 2 ) = Trp(B - TrpB) 2 = TrpB 2 - (TrpB) 2 (24a) 

be the squared variance of B formed with respect to p, and 

(((AB) 2 )) = Trp(B - TrpB) 2 = TrpB 2 - (TrpB) 2 (246) 

be the corresponding squared variance of B formed with respect to p. The first of these two 
squared variances fluctuates from unraveling to unraveling; taking its expectation over the 
stochastic process we have 

E[((AB) 2 )} =TrpB 2 - E[(TrpB) 2 } 

(25a) 

=Tr pB 2 - (TrpB) 2 + C , 

with C a correction term given by 

C =(TrE[p\B) 2 - E[(TrpB) 2 ] 

(256) 

= - E[(TrpB - TrE[p]B) 2 ] < . 
Hence we have obtained the inequality 

E[((AB) 2 )] < TrpB 2 - (TrpB) 2 = (((AB) 2 )) , (26) 

in other words, the squared variance formed from p gives an upper bound to the expectation 
of the squared variance formed from p. These results, and those of Sec. 3, can be used to 
calculate bounds on the expected variances ^[((AXi^) 2 )]. When the effects of the driving 
terms d(t), d(t) can be neglected (or at least remain bounded), we see, for example, that at 
large times we have from Eq. (10) 

E[((AX lt2 ) 2 )] < Trp(t)Xl 2 - (Trp(t)X lt2 ) 2 ~ 2 W 4 t , (27a) 
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giving a large time bound on the mean squared stochastic fluctuations of Xi$. Similarly, 
we find (using the discussion following Eq. (16d)) that when the effects of the driving terms 
can be neglected, the leading large time variance of N is bounded by 

E[({AN) 2 )] <Tr p{t)N 2 - (Trp(t)iV) 2 = Trp(t) (aVaa + a f a) - (Tr p(t)<f a) 2 

(276) 

~2 V 2 aH 2 - ( W 2 t) 2 = rfaH 2 . 
Thus the root mean square variance in N, and the expectation of N, have the same time 
rate of growth. 

7. Perturbation analysis for stochastic fluctuations 

We conclude our theoretical analysis by developing a formal perturbation theory for 
solving the evolution equation of Eq. (5a) for the pure state density matrix p. Let obey 
the evolution equation 

dpW = ~[H, p {0) ]dt, (28a) 

which holds when there are no stochastic terms, and let us expand the solution p of the 
corresponding stochastic equation as 

p = p(o) + v /r ? p( 1 /2) +r?/) (i) + ... . (286 ) 

Inserting this expansion into Eq. (5a), and equating like powers of i] on left and right, we 
get the following stochastic differential equations for p^ 1 / 2 ) and pfi\ 

d fiW) =- l -[H, pW]dt + a[p(°\ [p(°\ a + a^]}dW t , 

dpU = - l -[H,pW]dt - \a 2 [a + a\ [a + a\ p<®]]dt (28c) 

+<j([p(°\ [p^ 2 \a + at]] + [pW, [p(°), a + a%dW t . 
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The first step in solving these equations is to eliminate the time evolution associated with 
the Hamiltonian term %uxva by defining, for any operator B, an interaction picture operator 
B 1 given by 

so that in particular 

a 7 = ae-^ , a 7t = ate™ . (296) 



Then in interaction picture, Eqs. (28a) and (28c) become 

dp m = _ l [h i iP m ]dt y 

I i 

d piW2) = _ l [h i^ piW^ dt + (T [^(o) > [/0 /(o) > ae -iu,t + a \ e i^ dWt ; 
n 

dp 1 ^ = _ Uh^p^dt - \a 2 [ae~ wt + cte™, \ae~ wt + aV<**, p I{ ^\\dt 



(29c) 



+ ( 7([p / (°), [p^ 2 \ ae-™ + aV"*]] + [p^' 2 \ [p I{0 \ ae~ wt + aV"*]])^ . 
Here h 1 = h J (t) denotes the interaction picture form of the oscillator driving terms in the 
Hamiltonian, 

h T (t) = d{t)ah iu}t + d{t)ae- iujt . (29d) 

We can now deal with the h 1 term in the equations of motion by introducing an operator 
C/ / (t) that obeys the differential equation 

which, using the definitions of Eq. (16a), can be explicitly integrated to give 

U^t) =exp ^-^D(t)a^j exp (- l -D(t)a^j exp (J^ dud{u)e tu;u D{u)^ . (306) 
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We can now use U 1 as an integrating factor to integrate Eq. (29c), giving finally explicit 
formulas for p I( -°\ p 1 ^ 1 ^ 2 ^, and p I( ^\ 

p^(«)(t) =C/ / (t)p / W(0)c7 / (t)t , 

p%) = a f U I (t)U I ( S y[p m (s),[p I ^(s),ae- iu)S + (J e^U 1 (s)U 1 (tf dW s , 

JO 

p^it) =--a 2 / C/ 7 (t)C/ 7 (s)t[ae- ia;s + ate™, [ae~ iws + a^e iu)S , p m (s^U 1 (s)U 1 (t^ds 
2 Jo 

+a J V it) U 1 ( s) t ( [p /(0) (s) , [p /(1/2) (s) , ae"^ s + a f e iws ] ] 
+ [/3 / ( 1 / 2 )(s), [p / (°)(s),ae-^ + a^e^U 1 ^U 1 (tydW s . 

(31a) 

These equations give terms in the expansion in powers of ^Jr\ of the interaction picture 
quantity p 1 ; to transform back to the original Schrodinger picture, one uses the inverse of 
Eq. (29a), 

taking B to be successively p^ ' 1 / 2 - 1 ). From the results of this calculation, one can in prin- 
ciple compute the quantum expectations TipB corresponding to different unravelings of the 
stochastic process. 

We see from Eq. (31a) that the expression for p^ 1 / 2 ) involves a stochastic integration 
over dW s with a non-stochastic integrand, and so as expected we have E[p^^(t)] = 0. 
Additionally, we note that the expression for p^> contains two integrals, an ordinary integral 
involving an integration over ds, and a stochastic integral involving an integration over dW s . 
Since the integrand of the latter contains only stochastic quantities depending (through 
p J (V 2 )( s )) on dW t for t < s, the stochastic expectation of the integral over dW s is zero, and 
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so we have 

E[p I{l \t\) = -\a 2 f U\t)U\s)^[ae'^ s + a^e^\[ae- ws + ah w \p I ^\s)]}U\s)U\t) ] ds , 
2 Jo 

(32) 

which writing p 7 W(s) = Eft^^s)] gives the first term in the perturbation expansion for the 
ensemble expectation density matrix p(t) obeying Eq. (5b). 

Let us now use the results, £[p (1/2) (t)] = and E[p^(t)] = pW(t) to interpret the 
inequality derived in Sec. 6. Inserting the expansion for p into the definition of Eq. (24a), 
we have 

((A£?) 2 > = Tr(p(°) + v/??p (1/2) + rjp^)B 2 - (Tr(p^ + ^p^ + VP {1) )B) 2 . (33a) 

Taking now the expectation of this equation, we get for non-stochastic operators B, 

E[((AB) 2 )} = Tr(p(°) + rip^)B 2 - (Tr(p<® + r]P (1) )B) 2 - r]E[(Tip^ B) 2 ] + 0(rf ) . (336) 

But comparing now with Eq. (24b), we see that this is just 

£[<(A5) 2 >] = «(A5) 2 » - V E[(Tr(p^B) 2 ] + 0( V 2 ) , (33c) 

in agreement with the expansion of the inequality of Eq. (26) through terms of first order 
in T), and giving us insight into why the inequality takes this form. By writing Eq. (31a) for 
pi 1 / 2 ) in the form 

pW\t)= [ t p(s,t)dW s , (34a) 
Jo 

with P(s, t) denoting the integrand in Eq. (31a), and using the Ito isometry given in Appendix 
A, the stochastic expectation in the final term in Eq. (33c) can be explicitly evaluated as an 
ordinary integral, 

E[(Ti(p^B) 2 }= f ds(TiP(s,t)B) 2 . (346) 
Jo 
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8. Estimates for gravitational wave detection 
and nanomechanical oscillator experiments 

Let us now use the results of the preceding sections to make estimates for precision 
experiments involving monitoring of harmonically bound or free masses. We begin by col- 
lecting the relevant formulas. For the harmonic oscillator, we have seen in Eq. (9b) that the 
double expectation of the occupation number N = a^a has a secular growth given by 

({N)) ~ r]G 2 t. (35a) 

Since by our definition of Eq. (2a), a = (h/2muj)^ , we have [10] a = AX S q L , with AX S ql the 
so-called "standard quantum limit" for a conventional amplitude-and-phase measurement of 
Xi or X 2 , and so we can rewrite Eq. (35a) as 

((iV» ^(AX 5QL ) 2 t . (356) 

We have also seen in Eq. (27b) that the right-hand side of Eq. (35b) also gives at large times 
an upper bound to the root mean square variance in N, 

E[((AN) 2 )]^ <r)(AX SQL )H . (35c) 

For the quantum nondemolition variables X 12 , we have seen in Eq. (10) that the double 
expectation is not influenced by stochastic reduction or decoherence effects, 

6{{X lt2 )) = , (35d) 

while from Eq. (27a) we get at large times an upper bound to the root mean square variances 
in X h2 , 

£[((AX lj2 ) 2 >]5 < (27jf) V = (2r ] t^(AX SQL ) 2 ■ (36a) 

23 



This last equation has a similar form to the corresponding equation for a free particle, for 
which the standard quantum limit AqsQL in a position measurement is given [10] by 

Aq SQ L = (ht/m)* , (366) 

so that at large times we have from Eq. (13c), and the fact that STrpq = 0, 

£[((Ag) 2 )]^ < (oTrpg 2 )^ = ( V t/3)ht/m = (W#(Ag SQi ) 2 • (36c) 

These equations will form the basis for our analysis of experiments in which oscillating or 
free masses are monitored. Since we are making only order of magnitude estimates, we shall 
neglect numerical factors of order unity (such as the factor of 3 arising from generalizing 
from one to three dimensions), quoting all answers as powers of 10. 

To make estimates, we shall need values of both the stochasticity parameter r] and the 
elapsed time t. The value of 77 depends on the stochastic reduction model under consideration. 
In the GRW model [11] and also the QMUPL model [11], rj = rj N, with r] ~ lO^s^m" 2 
and with N the number of nucleons that are displaced in the measurement. For the CSL 
model, one has [8,11] 77 = ^ S 2 D 2 {a / -k)^ , with S the side length (for a cube of material), 
D the density and 7 ~ 10 -30 cm 3 s _1 and a ~ 10 10 cm -2 parameters of the model. We shall 
assume a nucleon density of D ~ 10 24 cm~ 3 , and shall ignore the geometry dependence of rj 
by eliminating S in terms of D and the nucleon number iV by writing S* 2 = (N/D)%, giving 

77 = 7 iV§ J D I (a/7r) 5 . (37a) 

For the elapsed time t we shall take the inverse of the noise bandwidth frequency F = 00/ (2Q), 
with Q the quality factor, for the nano mechanical resonator experiment, and the inverse of 
the low frequency limit of the sensitive range for the gravitational wave detector experiments. 
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Our reasoning here is that if accumulation of a small stochastic effect takes longer than the 
time estimated this way, the effects will be hard to distinguish from accumulated effects of 
the noise that sets the low frequency limit of the detector. 

The first experiment that we shall consider is the nanomechanical resonator reported 
by LaHaye et. al. [15], which uses a 19.7 MHz mechanical resonator containing ~ 10 12 
nucleons, corresponding to AX S ql ~ 1CT 14 m, and which has a noise bandwidth F = 903 Hz. 
For the GRW and QMUPL models, we have r] ~ 10 10 s _1 m -2 , giving an accumulated {{N)) 
in time F^ 1 of 10 -21 , and a root mean square expected deviation in X ij2 of ~ 10~ w AX S ql- 
In the CSL model, i] ~ 10 19 s -1 m -2 , giving an accumulated ((N)) in time F _1 of 10~ 12 , and 
a corresponding root mean square expected deviation in X 12 of ~ 10~ 6 AX S q L . 

The second experiment that we shall consider is the upgraded version of LIGO (the 
Advanced LIGO Interferometers), which monitors a quasi-free mass of 40 kg ~ 10 28 nucleons, 
and has a sensitive range extending down to F ~ 70 Hz. In the GRW model this apparatus 
has i] ~ 10 27 s _1 m -2 , while in the CSL model the value of this parameter is i] ~ 10 30 s _1 m -2 . 
The standard quantum limit of Eq. (36b) for position measurement over a time interval 
t = (70Hz)^ 1 is Aq SQL ~ 10~ 19 m, and we find that the root mean square stochastic 
deviation in the coordinate q over this time interval is bounded by ~ 10~ 7 AqsQL in the 
GRW model, and by ~ 10~ 5 Aq SQL in the CSL model. 

The third experiment that we consider is the projected space-based Laser Interfer- 
ometer Space Antenna (LISA) [17], which will monitor the positions of 2 kg masses to an 
accuracy of 10~ n m, and which will be sensitive to frequencies down to 10~ 4 Hz. From 
Eq. (36b), the standard quantum limit corresponding to a 2 kg mass and t ~ 10 4 s is 10~ 15 
m, in other words, this experiment will achieve a position accuracy of around 10 4 times 
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AqsQi- For the CSL model, the corresponding root mean square stochastic deviation in 
the coordinate will be of order lOOAqsQL, which is still a factor of 100 smaller than the 
observable displacement. 

We see that in nanomechanical oscillator and Advanced LIGO experiments, pre- 
dicted stochastic reduction effects are at least a factor of 10~ 5 below the relevant standard 
quantum limits, and so are presently far from being detectable. The situation is better 
for LISA, where the stochastic reduction effect is predicted to be two orders of magnitude 
larger than the standard quantum limit, but still two orders of magnitude below the design 
position sensitivity. Even though these experiments are not expected to observe an effect, 
they will place useful bounds on the stochasticity parameter r\. Trying to do better will be 
a challenging goal for future experiments; clearly, the key will be achieving a much larger 
accumulation time t, corresponding to a greatly reduced noise bandwidth F for the nanome- 
chanical resonator, or a greatly reduced lower frequency limit for the gravitational wave 
detectors. We note in closing that when Eq. (5b) is used as a model for environmentally 
induced (as opposed to postulated intrinsic) decoherence effects, the appropriate value of rj 
may be much larger than in the above estimates, and so in this case the effects for which we 
have obtained theoretical formulas may lie within reach of current experimental technique. 
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Appendix A: Basic Ito Calculus Formulas 

The stochastic differential dW t behaves heuristically as a random square root of dt, 
as expressed in the Ito calculus rules 

dW? = dt , dW t dt = dt 2 = . (Al) 

As a consequence of Eq. (Al), the Leibniz chain rule of the usual calculus is modified to 

d(AB) = (dA) B + AdB + dAdB . (A2) 

Applying these two formulas to the definition 

p(t) = \A)(A\ , (A3) 

the stochastic equation of motion of Eq. (4b) for \ip t ) is easily seen to imply the equation of 
motion of Eq. (5a) for p(t). Because the Ito differential dW t is statistically independent of 
the variables at and before time t, the final term in Eq. (5a) vanishes when the stochastic 
expectation E[p] is taken, giving Eq. (5b). Using this statistical independence of dW t , and 
Eqs. (A1)-(A3), we can get a useful formula for the expectation of a product of stochastic 
integrals. Consider the expectation 

fit) = E[ f dW u A(u) f dW u B(u)} , (A4) 
Jo Jo 
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which by Eq. (A2) has the differential 

df(t) = E[dW t A(t) [ dW u B(u) + [ dW u A(u)dW t B(t) + A(t)B(t)dt\ = E[A(t)B(t)\dt . 
Jo Jo 

(A5) 

Integrating back using the right hand side of Eq. (A5), we get 

E[ I dW u A(u) [ dW u B(u)} = [ duE[A{u)B{u)} , (A6) 
Jo Jo Jo 

a formula called the ltd isometry. 

Appendix B: Connection to the Lindblad Evolution Equation 

The most general completely positive density matrix evolution equation is given by 
the form studied by Lindblad [18] and Gorini, Kossakowski, and Sudarshan [18], generally 
referred to as the Lindblad equation, 

When Lj is self-adjoint, so that Lj = L], the summand in Eq. (Bl) reduces to the form 

which corresponds to the decoherence equation studied in the text when we take Lj = q, and 
more generally leads to a solvable oscillator model when Lj = Ciq + c 2 p (with self-adjointness 
requiring real c^.) These two cases correspond respectively to repeated environmental (or 
intrinsic, in the case of reduction models) measurements of the system coupling to q or to 
ciq + c 2 p- 

Dissipative equations in the Lindblad context are generated by taking Lj to be non- 
self-adjoint. For example, if we take Lj = a for the harmonic oscillator, we get additional 
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terms that cannot be represented as double commutators, as seen from the identity 
4<r 2 [apa t - ^{a i a,p}} = - ^[q, [q,p\] - — ^ [p, [p,p\] 

and a similar identity in which a is interchanged with a) and the i on the right hand side is 
replaced by —i. The paper of Salama and Gisin [14] studies a dissipative Lindblad equation 
with a term Lj oc a, while the papers of Isar, Sandulescu, and Scheid [14] and Karrlein and 
Grabert [14] study a non-Lindblad master equation with a single dissipative term propor- 
tional to [q, {p, p}]. 
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